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Abstract 

We study the competition between non-magnetic (dimer or flux) 
states and short-range antiferromagnetically ordered RVB states in the 
t-J model and present a finite temperature phase diagram. We show 
that, for a wide range of temperatures and dopings, the stable phase 
results to be a state in which both the flux and the RVB parameters 
are nonzero. 

PACS: 71.27, 74.20. 



In this paper we want to extend the analysis of the phase diagram of the 
t-J model on a square lattice, started by some of the authors in ref. 0, where 
only the so called non-magnetic states (such as the dimer and the flux phase) 
were considered. 



Our starting point is the t-J Hamiltonian, as originally proposed in 

J 

4 



n = t5 Y,J2 C L C i« - 4EE C la C J« C ]/3 C ^ , (i) 



where (ij) stands for a sum over ordered nearest-neighbor sites. This Hamilto- 
nian incorporates the below-half-filling constraint only in average: the hopping 
coefficient t of the first (kinetic) term gets renormalized by the doping factor 5 
fl. At half filling (5 = 0) the Hamiltonian reduces to the well-known Heisenberg 
Hamiltonian, with the spin operators written in a fermion representation. 

As candidates for the ground state of the t-J model, the RVB and the flux 
states have been mainly considered in literature. The former, first proposed by 
Anderson and coworkers 0, is characterized by a BCS-type order parameter 
of the kind A y - = (c^Cji — c^c^), which gives short-range antiferromagnetic 
correlations. The latter || corresponds to a nonmagnetic complex order pa- 
rameter Xij — (Ea=t,i c la c ja)> whose phase gives rise to a nonzero magnetic 
flux threading the elementary (square) plaquettes of the lattice in a staggered 
way. It is well known |3J] that the Heisenberg Hamiltonian, and hence the t-J 
model (H) at half filling, is invariant under local SU(2) transformations. One 
can exploit this symmetry H to show that the so called mixed-RVB and the 
flux phases are degenerate since they are related via one such gauge transfor- 
mation. Thus, these apparently different phases actually represent the same 
physical state. Away from half filling, this SU(2)-gaaige symmetry gets ex- 
plicitely broken by the kinetic term. This means that the RVB and the flux 
phases are no longer equivalent so that the question of a possible competition 
between these two phases arises. 

To start with, we want to reformulate the theory in a way that enables us 
to keep always track of the (broken) gauge symmetry. This can be done by 
introducing the 2x2 matrix operators: 



*< = 1' ^ (2) 

and by defining: 

$y = *i*t (3) 



§For a detailed discussion of how to derive this formula in the strong-coupling limit of 
the Hubbard model at or below half-filling see Q] and references therein. 
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Then, with some algebra, one can easily show that ([I]) can be rewritten as: 

H = f*5>[*s*«] - ^EM4 $ d > ( 4 ) 

(y) (ij) 

where C3 is the third Pauli matrix and tr stands for the trace over 2x2 
matrices. It is immediate to check that the second term in the Hamiltonian is 
gauge invariant under the following S , f/(2)-action: 

$ij — ► Ui^ijUj , (5) 

Qi being any element of SU(2) possibly site-dependent, while the hopping 
term explicitly breaks such symmetry. Both the RVB and the nonmagnetic 
phases of the t-J can be studied in this context. Indeed, in the saddle point 
approximation, the mean value of the fields $j/s is a function of Ay and Xiji 
namely: 

where * denotes complex conjugation. Bearing this in mind, we study the t-J 
Hamiltonian by representing the partition function in the grand-canonical en- 
semble by means of a functional integral over Grasmann fields representing the 
fermionic operators and by decoupling the four-fermion Heisenberg interaction 
via a Hubbard-Stratonovich transformation over the auxiliary SU (2)-fields 

U = ~ X * ij Aij I (7) 

Following [Q] and || we will assume invariance under translations along the 
lattice for the parameters A's and along the diagonal for the x' s - We are then 
left with a total of six independent parameters, A x<y and Xi, 2,3,4, as shown in 
Fig. |. 

Due to the assumed translational symmetry, it is convenient to introduce a 
different notation for the fermionic fields living on odd and even sites (empty 
and full dots in the figure), say ^'s and 0's. If we now go to momentum space 
- with k belonging to the Reduced Brillouin Zone (RBZ) and u n denoting the 
fermionic Matsubara frequencies - the partition function in the static approx- 
imation (i.e. the auxiliary fields being time independent) can be rewritten 
as: 



[»«<,,„ 2>VwJ P>««,» P*ba,„] / PXj v Xi\ \V A* T>Aj 
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Figure 1: 

The unit cell of the square lattice and the definition of 
(a) the RVB and (b) the flux order parameters. 
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where 



A fc = ( Xl+ t5)e^ + ( x ; + t5)e-^ + (x3 + ^)e-^ + (xl + ^)e^ (9) 
Afc = 2(A X cos fc x + A y cos (10) 

and M is the number of lattice sites. The effective action in the functional 
integral ([§]) is now quadratic in the Grassmann fields ip, <j), and it is actually a 
sum of independent quadratic forms acting separately in the four-dimensional 

subspaces spanned by {ip ka ,n, 4>ka,n, V-m -n> 0-k,a,-n)> for k in tlie RBZ > a =t> 1 
and c<j n = (2n + where they are represented by the matrices: 



( V 

Afc 



-A* 





-Afc 
-H 
-At 



-Afc \ 


-Afc 
-fi J 



V -Afc 

Integrating then over the Grasmann fields, we arrive at the result: 
Z = [[VxPxj][V&P^]exp{-m 



(12) 



n 



2M 



EN 2 + 2 E |a, 

3=1 



s, 



1 



l 



s, 



eff 



eff 



-ILOn 



(3 log E{) 



(13) 
(14) 



keRBZ j=l 

I 



= 4 £ Eiog(i+^; 

^ keRBZ j=l 

where (j = 1, 2, 3, 4) denote the eigenvalues of the matrix ([□]) multiplied by 
(3. Such matrix being 4x4, it would be possible to give an analytic formula for 
the eigenvalues E{ and their dependence on the parameters /i, A^, A&, which is 
however very complicated and will not be presented here for the general case. 

From (p~3|) one can then derive the self-consistency equation that must hold 
for the parameters x's and A's in the mean field approximation: 

dfl 

- - (,' = 1,2,3,4) (15) 



{j = x, y) 



(16) 



which have to be supplemented by the equation fixing the chemical potential: 

(17) 



We have numerically investigated the solutions to this set of equations that 
also minimize the grand potential Q, finding that all solutions belong to one 
of the following classes: 



1) Mixed-wave solutions : 

Xl=X2 = X3 



Xi = X 



A x = A , Ay 



e' T A . 



2) Dimer solutions : 
Xj real, 



Xi >> |X2 1 



1X3 | 



|X4| 



A* = A, 







(19) 



Under this assumptions it is possible to give a nice formula for the eigen- 
values E J k and hence for the grand potential. 
The effective action for the mixed-wave solution reads: 



S, 



eff 

A fc 



-- y 



keRBZ 



log cosh -— — + log cosh --— - + 2 log 2 



(20) 



(/x±|A fc |) 2 + |A fc |2 

2[(x + t$) cos + (x* + t$) COS ky) 
2 A (cos k x + e lT cos k y ) 
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For the dimer solution one finds: 




- [log(l + e^-^)) + log(l + e^+^)) 




(21) 



Ek 



|Afc| 

(Xi + t5Y k * + ( X * 2 + t5)e- ik » + ( X 3 + t5)e-^ + (xl + tS)e 



It is among these classes of solutions that one has to look for the absolute 
minimum of the free energy F = Q + yuM(l — 5) as the temperature T and the 
doping 5 vary. But, before presenting the phase diagram in the 5 — T plane, 
we need to understand some more properties of such solutions. 

1) The solution we call mixed- wave can actually represent different phases 
already studied in literature, depending on the relative values of the two 
independent complex parameters x an d A: 

a) Uniform phase: it corresponds to A = and x rea l- As we will 
see this is the stable phase at high doping and high temperature. 

b) Flux phase: it corresponds again to A = but now x nas a 
nonvanishing imaginary part. If we write x = \x\ e ^ we have that 
40 might be seen as a magnetic flux threading the elementary square 
plaquette. One can check that, at fixed T, the value of the phase 
increases while S decreases, going from = at high doping 
(uniform phase) to = 7r/4 at half filling (5 = 0). 

c) RVB-d phase: it is obtained when x — and r = it. The actual 
value of the parameter A depends, at fixed T, on the doping 5. 

d) RVB-mixed phase: it is given by x — and r = n/2. Again the 
value of the parameter A depends on the doping 5. 

e) It is important to notice that, at half filling, the 7r-flux and the RVB- 
mixed phases are equivalent via an SU (2) transformation, and hence 
are degenarate in energy and actually represent the same physical 
state. It is only away from half filling that these two phases become 
independent. Still, one expects them to remain very close in energy, 
at least for low doping. 

2) In the dimer solution, the relative values of the two independent order 
parameters Xi an d X2 change with 5, for fixed T. In particular, at any 
temperature, X2 = at half filling, so that only one link variable per 
plaquette is nonzero, giving rise to the so called staggered dimer config- 
uration. 

As for solutions of type (1), our analysis has shown that none of the above 
mentioned "pure solutions" is the actual absolute minimum of the free energy. 
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At not too small values of 5, the latter is indeed minimized by a phase in which 
both the x an d the A parameters are different from zero. Precisely we obtain 
a sort of mixture between the uniform and the RVB-d phases, which we would 
like to call the mixed phase^ 

X real and r = n . (22) 

The values of \x\ and |A| change with T and 5. Above a critical temperature 
T c = J/8 we find A = and hence recover the uniform solution. At any 
T < T c , there is a critical doping S c above which the stable solution is again 
the uniform phase. Below S c , the parameter A increases continuously (second 
order transition) while S decreases, in such a way that for 5 = one also 
has \x\ = |A|. It is interesting to notice that the solution we find at half 
filling, where the SU(2) symmetry is restored, is indeed equivalent via a gauge 
transformation to the RVB-mixed phase, or to the 7r-flux phase. This means 
that the introduction of the kinetic term breaks the S{7(2)-degeneracy in favour 
of a state where all the x' s are equal and real, while the A's differ just for a 
phase of it. 

At low dopings, it is the dimer phase that becomes stable. The transition 
from the dimer to the mixed phase is a first order transition. 

The phase diagram is shown in Fig. |2] for t/J = 1. For different values 
oit/J one obtains diagrams that look qualitatively the same: what changes 
slightly is only the position of the transition lines. 

We have also investigated the stability of the mean field solutions around 
the lines of transition, by looking at the behaviour of the chemical potential 
/j as a function of doping S, since an instability towards phase separation is 
naled by |f > (see the first paper in M for a discussion). 

While no instability develops at the mixed/uniform line of transition, we 
find that the coexistence of the dimer and the mixed phases is favoured at all 
temperatures < T < T c around the original line of transition. Such effect 
is pretty strong especially at low temperatures, so that the pure dimer phase 
survives only very close to half filling. This is shown in Fig. || 

We would like now to compare these results with our previous work M. 
Even if the starting Hamiltonian and the formalism we adopted are the same, 
and despite the fact that in the present paper we have enlarged the class of 
allowed mean field solutions, we find here a much simpler phase diagram. In 
particular we notice that the flux phase is nolonger present, at any value of 



the doping and temperature, being substituted by the mixed phase (|22|). This 
is in agreement to and extends the findings of Zhang at zero temperature 0. 
Also, our mean field Hamiltonian and self-consistency equations coincide, but 

**This phase is simply referred to as the "dimer phase" in |7|, ^ 
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0.14 




Doping (5) 

Figure 2: 

The mean- field phase diagram of the t-J model. 



for some normalization coefficients, with the one given in |§ and obtained 
through an auxiliary-boson approach, originally formulated by Kotliar and 
Ruckenstein |J for the Hubbard model, in which an empty site is described 
by a bosonic ceation operator b\, whose density gets to be proportional to the 
doping in the meand field approximation: (b\bj) = 50. 

On the contrary, the authors of ref. || find some wide region where the 
flux phase is stable and a small area of the phase diagram where there is a 
coexistence of the flux and the superconducting phases. They represent the t-J 
model via the so-called Hubbard X-operators, which allow for the enforcemeent 
of the below-half filling constraint, and perform a 1/N espansion to extrapolate 
then to the case N = 2. 

Our phase diagram is also substantially different from the one presented by 



Ubbens and Lee in [10] and obtained within a slave boson approach that differs 



from the one mentioned above || by the way the hopping term is treated. The 



free energy and hence the self-consistency equations found in [10| have an 
additional contribution coming from all the (infinite) bosonic states. As a 
result, the mean- field phase diagram shows a finite region, close to half filling, 
where the flux solution (either tt or staggered) is favoured against all other 
phases. A similar result is obtained, at zero temperature, also by Sheng et al. 
in ||11|. They work again within the slave boson approach by introducing now 



* Parenthetically, we notice that this explains also why the hopping Hamiltonian, which 
is written as t J2 a c ia c ja^jbi in this approach, becomes t5 J2(ij) S Q c la c j°< m ^ ne mean 
field approximation as we wrote in dl]). 
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Figure 3: 

The phase diagram of the t- J model, including the 
occurrence of phase separation. 



an additional statistical gauge field that describes the hard-core nature of holes. 
In this case the free energy gets minimized, at low dopings, by a staggered flux 
phase combined with a nonzero statistical field yielding a uniform flux per 
plaquette. 

As we have also noticed in the second paper of ref. our calculations (a la 
Anderson) and the slave boson techniques seem to lead to very different results, 
in particular close to half filling. We are investigating why this discrepancy 
occurs. Even if we do not have a final answer to this question, we believe 
that the problem has to do with the way the below-half filling constraint is 
treated within the different approaches. Indeed, at the mean field level such 
constraint is respected only in some average sense and differences in the way 
this is implemented can result in different phase diagrams, especially close 
to half filling when the high density of electrons enhances the probability of 
including configurations with double occupied sites. 

Finally, we want to comment on the possibility that, at a fixed temperature, 
phase separation might disappear when J/t is lowered below a critical value 
J c /t. In the second paper of ref. ||J, we studied how this happens when 
considering non- magnetic phases only, showing that J c /t varies linearly with 
T according to J c / 1 « 8T. This results is confirmed in our present calculations, 
where we have also included RVB states. There is a difference however: now 
phase separation occurs as soon as the temperature is lower than the critical 
temperature T c = J/8, whereas in our previous paper at temperatures close to 
the critical one only pure phases were present. 
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